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SUMMARY 

The feasibility of a laser-driven MHD generator, as a candidate receiver 
for a space-based laser power transmission system, was investigated. 

An extensive literature search of research on MHD generators and laser- 
produced plasmas was carried out. The MHD generators were tabulated accord- 
ing to characteristics such as the energy source, working fluid, generator 
type, flow rate, temperature, electrical conductivity, power density, gener- 
ator dimension, efficiency, magnetic field strength, seed material, type of 
cycle, and operating mode. Laser-produced plasma and laser plasma inter- 
actions were tabulated with respect to plasma temperature, laser type and 
energy, plasma conductivity, absorption of laser radiation, flow velocity, 
carrier gas, and seed material. 

On the basis of reasonable parameters obtained in the literature search, 
a model of the laser-driven MHD generator was developed with the assumptions 
of a steady, turbulent, two-dimensional flow. The assumptions used in this 
study were based on the continuous and steady generation of plasmas by the 
exposure of the continuous wave laser beam thus inducing a steady back pres- 
sure that enables the medium to flow steadily. The model considered here 
took the turbulent nature of plasmas into account in the two-dimensional 
geometry of the generator. For these conditions with the plasma parameters 
defining the thermal conductivity, viscosity, electrical conductivity for 
the plasma flow, a generator efficiency of 53.5 percent was calculated. If 
turbulent effects and nonequilibrium ionization are taken into account, the 
efficiency is 43.2 percent. 

The study shows that the laser-driven MHD system has potential as a 
laser power receiver for space applications because of its high energy 
conversion efficiency, high energy density and relatively simple mechanism 
as compared to other energy conversion cycles. 
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INTRODUCTION 


The advantages of using a laser to transmit power in space is based on 
three features of the laser: 

1. The laser beam can be transmitted over long distances 
without appreciable attenuation or divergence, 

2. The laser provides a high source intensity at the 
receiver, 

3. The laser does not require physical contact between 
energy source and power generator. 


While the laser appears to be an advantageous means of transmitting en- 
ergy in space, the means of beam generation and beam conversion to a more use- 
ful form of energy (i.e., electricity or propulsion) are not well defined. 
Converter systems, in particular, must meet the requirements of high con- 
version efficiency at a high energy density while remaining small, light 
weight and simple. One converter system which may meet the requirements 
is the laser-driven MHD generator which is shown in Figure 1. The plasma 
is produced by focusing the laser beam into the plasma production chamber. 

At sufficiently high intensity, breakdown will occur in the gas medium 
producing a plasma. Once the plasma is established, it will absorb the 
laser radiation which will heat the plasma. Although the laser driven MHD 
generator is an efficient candidate energy conversion system for space ab- 
plication, the absorption of transmitted beam energy by the participating 
medium in the plasma production chamber becomes the key factor in deter- 
mining the generator efficiency. With the proper plasma conditions such 
as laser peak power, gas pressure, gas species, focal volume density, 
and plasma temperature, absorption of the laser beam can reach 80% with 
Nd long pulse laser light of approximately 10^® W/m^ and 65% with CO 2 
long pulse laser light of roughly 10^ W/m^ (Ref. 1). Figure 2 shows 
the estimated efficiencies of a laser-driven MHD generator subsystem. 


This preliminary study of the laser driven MHD generator includes a 
literature survey and a development of a simplified theoretical model for 
the system. The literature survey was conducted to establish realistic para- 
meters, such as temperature and density, for laser produced plasmas and also 
to determine those design parameters of MHD channels which are affected by 
the plasma conditions such as plasma temperature and density. Both plasma 
and liquid metal MHD generators were included in the literature survey. 


The primary purpose of the study was to identify certain characteristics 
of existing MHD generator systems so that these characteristics could be used 
in the development of a prototype design for a laser driven MHD generator. 
System characteristics such as the energy source, generator Lype, working 
fluid, working fluid flow, temperature, electrical conductivity, power den- 
sity, generator dimension, magnetic field, cycle features, seed material and 
mode of operation were considered. 


1. Plasma MHD 


The characteristics of plasma MHD generator systems are given in Table 
I. The systems can be divided into four groups: 1) shock-driven, 2) arc, 

3) combustion (including coal burning), and 4) explosive-driven (including 
rocket and detonation). The primary difference between these systems and 
the laser-driven MHD generator is the method of plasma production, although 
some differences in the flow of the plasma through the channel occurs. The 
shock or explosive driven MHD generators, for example, are characterized by 
quasi-adiabatic wave propagation through the channel and a rapid decay of 
the plasma after the wave has passed. In the laser-driven MHD generator, on 
the other hand, the plasma flow will be exposed to the laser radiation (so 
long as the critical charge density which produces optical reflection is not 
reached), and heating of the plasma throughout the channel will occur. The 
flow in the laser-driven MHD channel is expected to be an unsteady, turbu- 
lent flow with a heat source. The power density in a MHD generator depends 
on the average channel velocity, the magnetic field intensity, and the elec- 
trical conductivity. 


2. Liquid Metal MHD 

The characteristics of liquid metal MHD systems are given in Table II. 

The conceptual and experimental works done so far are directly applicable 
to the laser-driven liquid metal MHD system, since the configuration of the 
conventional liquid metal MHD systems is the same except for the laser beam 
receiver. The laser energy is supplied through the optical system in which 
the beam energy is converted into the driving power of the liquid metal 
(two-phase flow) . 

Liquid metal MHD systems usually deal with low temperature energy sources. 
Hence, they do not have the hardware related problems, such as melting, that 
may take place in the plasma MHD generator. The corrosion and the separation 
of gas from liquid (in the case of two-phase flow) are the main problems to 
be solved. 

Table III gives the alkali metal parameters such as the ionization 
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potential the percent ionization, and the absorption length at 2500K and 
10333 Pa(l atm). 

3. Laser Plasma Interaction 


The plasma is produced by a gas breakdown after threshold irradiance has 
been achieved. This gas breakdown process proceeds in the following steps: 

(1) the production of the initial ionization, and (2) the subsequent cascade 
by which the ioniztion grows and the shock wave propagates, and dissipates. 

After the initial ionization is produced, its growth becomes the domi- 
nant process. Following a small amount of ionization, free electrons absorb 
photon energy by inverse bremsstrahlung. When an electron has gained enough 
energy, it can ionize an additional atom in a collision. The electron is 
then replaced by two electrons with lower energy in the free electron con- 
tinuum. Both electrons then absorb energy by inverse bremsstrahlung, and 
cascading of the ionization occurs. This cascade process, fed by the ab- 
sorption of laser light in the inverse bremsstrahlung process, is the mech- 
anism which produces the growth of the ionization. 

Breakdown threshold of the gas, as a function of intensity, depends on 
the focal volume. As the focal volume becomes smaller, losses, either by 
diffusion of the electrons out of the focal region or by radiation, limit 
the build up and increase the threshold. The cascade process proceeds more 
rapidly, for a given irradiance, with a larger focal volume. To maintain 
the growth of the cascade for a stable plasma within the focal volume, the 
following criteria should be met: 

(1) large focal volume 

18 2 

(2) high beam flux density (over 7 x 10 W/m ) 

(3) continuous wave beam flux (CW laser) 

(4) high gas density 

(5) low boundary effects. 


The power output requirements of the laser as an energy source for a 
laser-driven MHD generator can be chosen by the focal volume of the plasma 
chamber, gas pressure, gas species, MHD channel geometry, etc. 

Breakdown characteristics are shown in Figures 2, 4, 5 and 6. The gas 
breakdown threshold is a function of the gas pressure and focal volume. The 
breakdown threshold exhibits a linear relationship with the peak irradiance, 
but this relationship also depends on the species and pressures of the gas. 
For argon gas at pressures of 1 atm, 2 atm, and 38.2 atm (lines C, D, and 
E shown in Fig. 5), the breakdown has a different behavior at 1 atm than 
at 3 atm. At constant pressure of a participating gas, the absorption 
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of laser beam energy depends on the peak input power. Figure 6 shows the as- 
pect of beam transmission through the air at a pressure of 746 torr when a 
Ruby laser pulse is focused through a 0.0206 m focal length lens. Breakdown 
is also a function of the incident radiation wavelength and the density of 
the medium (see Fig. 7). The absorption of the laser beam occurs at high 
pressure (or density), and the time required for breakdown threshold must 
be short due to fast energy accumulations. At this state after breakdown, 
the dominant energy transfer mechanism (on a microscopic time scale) is 
inverse bremsstrahlung rather than conduction or diffusion. 

In Table IV, the relevant laser-plasma interaction parameters obtained 
from the literature survey are tabulated in order of the laser sources, tar- 
get medium, breakdown laser power, plasma temperature, plasma electron den- 
sity, dimension of the plasma, and propagation velocity. 

The laser considered in the tabulation of Table IV are the 10.6 ym 
CO 2 , Ruby pulse, and Nd-glass lasers. The media were air, argon, helium, 
nitrogen, deuterium, and plastic pellet at various pressures. 
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The MHD system selected here for study is a plasma MHD generator. The 
model of a selected MHD generator is set-up and simplified with the assump- 
tions that the system is under operation. 


1. MHD Channel Model 

The channel schematic selected for study is depicted in Figures 1 and 8 
and consists of a hydrodynamically, thermally stable turbulent flow of an elec- 
trically conducting and radiating gas, between electrode plates with a uniform, 
constant magnetic field applied in the positive x-direction. The channel side 
walls and the electrode walls have either a constant heat flux or a constant 
temperature. The physical properties of the gas are constant and the gas is 
in local thermodynamic equilibrium. The gas has a refractive index of unity, 
and scattering effects are negligible. 

•For the calculations, the geometrical dimensions and the boundary con- 
ditions for the channel must be considered. Table V shows the parameters 
which are necessary for calculating the MHD generator performance. The ini- 
tial and boundary conditions for the governing equations are determined from 
solutions to the equations for the laser-plasma interaction. For simulation 
purposes, the initial and boundary conditions can be obtained from experi- 
mental data. 

2. Energy Balance 

The thermal energy equation in this study considers the fully devel- 
oped, turbulent flow of an electrically conducting and radiating gas in a 
rectangular duct, where viscous dissipation and Joule heating effects are 
considered; but axial components of conduction and radiation are neglected 
(since their contributions to the velocity and temperature fields are small 
compared to those of the other terms in the equation). Thus the energy 
equation may be written in the following form: 
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The first two terms on the right hand side (RRS)‘ represent the net ther- 
mal energy transport due to molecular flow and turbulent transport with 
denoting the turbulent diffusivity of heat. The rest of the terms on the RHS 
are (in order shown) : the molecular and turbulent viscous dissipation with 

v t being the turbulent viscous dissipation of momentum, Joule heating with 
o e ff denoting the gas electrical conductivity, and divergence of the radiative 
flux. The turbulent transport quantities are assumed to vary across the chan- 
nel. The boundary conditions for the above equations are described by two 
possible cases: 


1) Constant wall temperature 

T * (y— 0) = T'(y- ± H/2) = T' 

' 1 

t ’ ( x=o > = t’ ( x = ± = r 



(2a) 



2) Constant heat flux 
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where q" and q" are the heat fluxes through the boundary. 

V w 2 .... ''' ■' • 

The initial conditions depend on the entry plasma characteristics. 
These can be defined by experimental measurement and treated as input vari- 
ables. Initially these conditions, listed in the bottom of Table II, were 
selected from the literature for the study. The boundary conditions in the 
z-direction totally depend on the initial conditions and on the time depen- 
dence of temperature at every nodal point. 

The following dimensionless quantities are introduced 


x = %, y = z - e 


L 


T' 

w 


T' T 
or 


V = — 2- 

t ■ * r u 


t = 


L/U 


Pr = 


Pr = — 

t a. 


Ec - 


c T ’ 

p w 


Ha = (-%^ j.- -4- 

;-.i, 2 m p o UB. 


e~ m 


N 


xk 
3 ’ 

4aT 

w 


T = xH, 

OX 


T = xW, Q 

OX X 


4t oT 
ox wx 


6 


4t 



oy wy 


where Pr and Pr t = molecular and turbulent Prandtl numbers, respectively. 

H 

Ha = the Hartmann number based on the half-channel width — , 

N = the ratio of molecular thermal conduction to radiation for 

a gas with an absorption coefficient k and optical thickness 

T . 

O 

Before nondimensionalizing the energy equation, it can be simplified by 
L'z >> an< * ^y ( s i nce the velocities perpendicular to the axial direction 
are much smaller in low Prandtl number fluids as compared to the axial flow, 
even in the boundary). Therefore, the 2nd and 3rd terms on the left hand 
side can be neglected. Applying the nondimens ional variables we have the 
energy equation 
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For steady plasma flow between the two parallel electrode plates, y-direction 
variaton only, no viscous dissipation, and the plasma with a dominant con- 
ductive loss in the channel, the above equation can be written in the simple 

form 
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The above equation, coupled with the following equations, is numerically 
solved by the finite difference method. The purpose of the computation 
was to show the temperature and the electrical conductivity profiles be- 
tween two electrodes along the MHD channel, to give the power output 
assuming a load to be applied, and to determine the generator efficiency. 

The power output density, considering the effective electrical conductivity 
is given by: 

2 

P = ~~ a U 2 B 2 K(1 - K) (5) 

i+s 2 eff 
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where 


K = ( -\/ 1 + 6 2 - 1) 3 2 


The current in the MHD channel, is 


J = -v/o « P 

Ay eff 


( 6 ) 


and the efficiency is 


n = 


B (1-K) 
K _1 + g 2 


(7) 


Velocity 

Kruger and Sonju (Ref. 2), employing the Karman-Pollhausen technique, 
estimated the wall shear stress and boundary-layer thickness corresponding 
to the semi-empirical velocity correlations proposed by Harris (Ref. 3). 

The local velocity normalized with the centerline value is evaluated from 

= (-y)* 5 |^6.154 + 2.457 In (Re t n) + F^—- n)J (8) 

where n is equivalent to y. 

Graphical results for the asymptotic friction coefficient presented by 
Kruger and Sonju (Ref. 1) may be approximated by 


= [lO. 
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where B = Ha/Re is the interaction parameter, 
is defined bv 


C f !< P U c L C f !, 
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where 0 < q < L. 


(B)J x 10 3 (9) 

The turbulent Reynolds number 

( 10 ) 


The function F^(^) is presented graphically by Harris (Ref. 3) and may 
be approximated by the following expression (Ref. 2): 
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for ; < 0.6 

F (;) = 2.502 + 21.930; - (6.359 + 53.747; + 649.535 c})' 2 (11a) 

and 


for ; > 0.6 

F (;) =-2.07 - 2.457 In (;) . (lib) 

1 2 

Ha 

where the parameter ; is equivalent to n- 


Near the wall, when n is small, velocities evaluated from Equation (8) become 
negative as a result of the logarithmic term. Therefore, in a manner similar 
to that employed in ordinary hydrodynamic (0HD) turbulent flows for the 
laminar sublayer, velocities are calculated utilizing the product of Re t *n up 
to the position where this product equals Equation (8). For the sublayer 
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where r^ is the hydraulic radius. For the rectangular cross section area of 
the channel, the hydraulic radius is 


_ LH 
r h 4 (L+H) 


(14) 


Therefore, the velocity in the sublayer is 


~ - 0.01757 ) 

c 



n- 


(13) 


Turbulent Viscosity 

For the viscous boundary and wall heat loss, the viscosity and thermal 
conductivity must be considered. In most practical MHI) applications, the 
flow is turbulent so that the transport processes are dominated by turbu- 
lent flow. Expressions for the turbulent viscosity for MHD flows are gen- 
erally based on those for turbulent viscosity for the MHD flow with modi- 
fications to account for such factors as the damping of the turbulent 
viscosity as the magnetic field is increased. 
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The OHD turbulent viscosity model used by Van Driest (Ref. 4) is modified 
by the Mei and Squire channel factor (Ref. 5), and is utilized with a multi- 
plicative magnetic damping function by Fiveland (Ref. 6). With these cor- 
rections , 


— = °‘ 5 - j i [l + 4K 2 Z 2 (1 - e Z/A ) 
v 1 + b n ( *- 

and the turbulent distance, Z, is defined by 

Z = nRe t . ( 17 > 



The magnetic damping function used by Fiveland (Ref. 6) is 


D 


-700 Ha 2 /Re 2 
e t 


(18) 


Thermal Conductivity 


In a plasma, unlike the case of viscosity, the internal structure of the 
colliding particles play an important role in determining the thermal con- 
ductibity. This is due to the fact that energy may be stored in internal de- 
grees of freedom such as rotation, vibration, and electronic excitation. In 
a mixture, which is in thermal equilibrium, particles recombine when they 
move against a temperature gradient and then release the energies of dis- 
sociation or ionization. However, the kinetic theory provides the simplest 
methods for estimating the transfer coefficients for a single component, mono- 
atomic gas. Denoting the particle concentration by n e , the mass of the par- 
ticles by m e , a nd the effective collision cross section for a solid sphere 
molecular model Q(T) , the coefficient of thermal conductivity for ionized gas 
is described by 


k = 


25 

16 


/ itnipicT 
Q(T) 



(19) 


The equation describes the thermal conductivity which would be valid if 
the composition of dissociated gas were frozen. That is, the processes are 
ideally faster than any chemical kinetic process. 

In Equation 19, the collision cross2sect ions , Q(T) , for rigid spheri- 
cal molecules of diam eter D is equal to nD 7 . This relationship can be 
used to estimate Q(T) for collisions between like molecules. 
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Electrical Conductivity 


The current density in the MHD channel is proportional to the electrical 
conductivity of the working medium. The formulation of current density, j , 
in the channel is given 

j = oU z B(1 - K) (20) 

where K is the generator coefficient and a is the electrical conductivity. 

The electrical conductivity is related to the electron density, n e , and the 
mobility, y, of the gas by a general form (Ref. 7): 

a = n ey (21) 


For MHD generators using alkali seeded noble gases as working fluids, non- 
equilibrium ionization occurs when Joule heating of the gas by the current 
causes the electron temperature, T e , to be higher than the gas temperature, 
Tg. To compute the electrical conductivity, taking into account the non- 
equilibrium ionization in MHD generators, the effective electrical con- 
ductivity rather than the scalar conductivity has to be used in the basic 
MHD equations. 

The effective conductivity is given by Zampaglione (Ref. 8) as 


a 


eff 


af(B - O 2 + (8P 2 1 

b[b + l (b 2 - i)] 


( 22 ) 


where B is the Hall parameter and £ is a plasma turbulence factor ranging 
from 0.5 to 1.0. The Hall parameter is 

B - uB . (23) 

The electron mobility, y, is given by 


V2 


(24) 


where V 2 is the collision frequency of electrons and neutrals. 


To calculate the electron density, n e , Saha's equation is modified because 
T e > Tg and the effective ionization potential of the seed is lowered by the 
Debye cloud. The modified Saha's equation for n e (Refs. 9, 10 and 11) is 



( 25 ) 
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K 1 = Z°^ ~ ^ " 3 - GX P [- 6 (V o " ^ )/<T e] • 

In these equations, F is the mole fraction of the seed, p is the total pres- 
sure, Z + is the electronic partition function of the seed ion, Zg is the 
partition function of the seed neutrals, h' is Plank's constant, V Q is the 
ionization potential, and y is the lowering of the ionization potential by 
the Debye cloud. The lowering factor of the ionization potential by the 
Debye cloud is defined by 


Y 


Z e 


4tt c p. 
o d 


(26) 


The Debye radius is given by 


p d = 


V 


e kT 
o e 

2 e 2 n 


(27) 


where z = 1 for atoms, 2 for + ions, and 3 for ++ ions, and e G is the permi- 
tivity of free space. 

The degree of ionization is given by 


a 


n c 


<n B + "Cs 1 


(28) 


where ng is the density of He or Ar. The collision frequency for He, as a 
function of a, is then defined by the approximation (Ref. 12), 

v 2 = [3.10 + (8085 a - 0.2264) 0,71 ] x 10~ 14 n He , (29) 

if (8085 a - 0.2264) > 0. Equation (29) becomes 
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if (8085 a - 0.2264) < 0. 


_i /. 

v 2 = 3.10 x 10 " n He* 


/ o r\ \ 




For Ar 

v 2 = ^0.53 + 0.641 x (10 4 a)°‘ 72 j x lO -14 n Ar 

3 

where n and n. are here in units of particle/cm . 
He Ar 


(31) 


The initial electron temperature at the entrance of the MHD channel can 
be defined by solving a set of equations describing the laser-plasma inter- 
action. The electron temperature in the MHD channel can be replaced by the 
plasma gas temperature which is implicitly computed by a set of equations 
for the MHD channel. In this simplified model, the entrance electron temp- 
erature was set at 2500 K. 


13 


RESULTS 


The breakdown threshold for plasma production by laser radiation depends 
on the medium pressure (Fig. 3), the focal volume (Fig. 4), the peak ir- 
radiance (Figs. 5 and 6), and the absorption band (Fig. 7). The growth of 
plasma beyond the breakdown threshold, however, depends entirely on absorption 
which falls into two different catagories, namely, short and long pulses 
(Ref. 1). In a single short pulse region, inverse bremsstrahlung comprises 
about 40% absorption. On the other hand, in a long pulse, absorption of over 
80% occurs up to the point where the transition from classical to anomalous 
behavior begins. Beyond the transition point, the absorption rapidly drops 
to about 50% due to the ion-acoustic turbulence, the Brillouin backscatter, 
the specular reflection, and the non-linear behavior of the plasma. The ab- 
sorption of high laser irradiance in the long pulse mode by an expanding 
plasma ball can be improved by a well-designed interaction chamber. The 
specular reflection and the scattering at the critical density surface of the 
plasma comprise about 40% of the total irradiance in the long pulse mode. 

Such losses are sensitive to the laser power level and can be partly re- 
trieved by considering the design of the optics and geometrical configuration, 
magnetically controlled plasma boundary, and by fabricating the plasma chamber 
with highly reflective surfaces. That is, well-designed optics and a geometry 
configuration with highly reflective surfaces can refocus the beam reflected 
from the plasma surface back to the plasma. The damage to the reflective 
chamber walls can also be alleviated by magnetically controlling the plasma 
boundary. In this case, the total absorption would be above 80% at the 
peak irradiance of a long pulsed laser. This is the optimized value of the 
conversion efficiency of the laser incident beam energy as shown in Figure 
2 . 


The characteristics of this absorption mechanism, as well as the crucial 
effect of the cold boundary, have to be accounted for in the plasma expansion, 
since the possibility exists that the plasma would decay before passing through 
the MHD generator. Applying a magnetic field to pinch the plasma radially un- 
til the plasma passes through the generator may solve the cold boundary prob- 
lem, or, on the other hand, alleviates the damage to the wall by the highly ex- 
panding plasma. The magnetic field may be used to thermally insulate the plasma 
from the cold wall. Since the cold wall is in contact with the plasma, it in- 
duces intabilities which enhance thermal conduction losses (Ref. 13). Table 
VI describes the laser sources, the medium, the breakdown threshold, the plasma 
temperature, the focal volume, and more. More research is needed to establish 
the dynamics and characteristics of the laser-induced plasma with respect to 
the irradiant power, the time to breakdown, growth, decay, the total absorption, 
and the medium pressure. How effectively the laser energy can be used to pro- 
duce a plasma in a chamber of improved geometrical features with specific 
radiative properties should also be investigated. 


The power density and the system efficiency of the MHD cycle is much 
higher compared to those of the conventional thermal cycles. However, the 
main problems in the development of MHD generators are the high temperatures 
required and the corrosive gas medium which can easily damage the electrodes 
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and wall of the generator. 

Though well-known concepts and well-developed analysis exist, the MHD 
generators still require intensive research for use in space. The existing 
analysis related to the conventional MHD generator may not be sufficient 
for the direct application to the laser-driven plasma MHD system because 
the energy source is in the form of a light beam. The light beam not only 
penetrates the plasma production chamber, but also extends the energy input 
through the generator, and thus may deposit the input energy unevenly 
throughout the plasma. 

From a simplified model computation, the temperature, velocity, and 
electrical conductivity distribution were obtained. These distribution 
curves from the center of the channel to the electrode at different points 
along the axial direction are shown in the Figures 9, 10, 11, and 12. 

The variables were non-dimensionalized so that the numbers in the figures 
are described with the scaled values. The velocity profile in Figure 9 
was obtained with the assumption of a turbulent channel flow. This ap- 
proach is reasonable, since the MHD flow usually has a turbulent flow 
profile. Based upon the velocity profile, the temperature and electrical 
conductivity for the medium were calculated. The temperature profiles in 
Figure 10 are scaled by the uniform wall temperature. By considering the 
electron mobility and plasma turbulences with the above velocity and temp- 
erature profiles, the effective electrical conductivity was computed and 
is shown in Figure 12. The electrical conductivity is significantly af- 
fected (as much as 55 percent) by the electron mobility and the plasma tur- 
bulence (which is also a function of temperature) . From the above computa- 
tional results, the average values of velocity, temperature, and effective 
electrical conductivity are obtained. These average values can be used to 
calculate the Hall parameter and the generator coefficient which are, in 
turn, used to calculate the power output density and the efficiency. The 
power output density and the efficiency considering the electron mobility 
and plasma turbulence was 0.2254 W/cm^ and 43 percent. These values can be 
improved by optimizing the design of the channel. For the same power out- 
put density, an efficiency of 54 percent is calculated based on the electron 
density of 1.315 x 10l9/ cin 3, the medium average temperature, 2500 K, and ig- 
noring turbulent effects and non-equilibrium ionization. Table VII shows 
che parameters of a MHD generator used for the two cases. 
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CONCLUSIONS 


In conclusion, the literature survey and the simplified model calculations 
show that based on its efficiency and power density, the laser-driven MHD sys- 
tem is practical and may have applications for future spacecraft energy system 
in space. 

From the simplified model the generator efficiency is 53.5 percent, if tur- 
bulence and non-equilibrium ionization are ignored. If these effects are taken 
into account, the efficiency is reduced to 43.2 percent. 

Further detailed theoretical and experimental analysis of the laser- 
driven MHD generator are necessary to realize a high potential for space ap- 
plication. 
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FIGURE 2. LASER-DRIVEN MHD SYSTEM EFFICIENCY BLOCK DIAGRAM 
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FIGURE 3. A COMPILATION OF THE EXPERIMENTAL RESULTS ON 
BREAKDOWN THRESHOLD AS A FUNCTION OF PRESSURE 
FOR A NUMBER OF GASES. (REF. 32) 
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FIGURE 5. BREAKDOWN TIME AS A FUNCTION OF PEAK IRRADIANCE 
AND PEAK ELECTRIC FIELD FOR A Q-SWITCHED RUBY 
LASER PULSE FOCUSED IN VARIOUS GASES. (REF. 54) 
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FIGURE 6. OPTICAL TRANSMISSIVITY OF AIR AT A PRESSURE OF 746 TORR AS A 
FUNCTION OF PEAK POWER IN A RUBY LASER PULSE FOCUSED BY A 
2.06 - cm FOCAL LENGTH LENS. (REF. 55) 




WAVELENGTH IA1 


FIGURE 7. BREAKDOWN THRESHOLD AS A FUNCTION OF 
WAVELENGTH OF INPUT RADIATION FOR Ar 
AT FOUR SELECTED PRESSURES. (REF. 56) 
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FIGURE 8. MHD CHANNEL 

J is the current density along the y direction 
B is the magnetic flux along the x direction 
U is the velocity along the z direction 
H is the height of the generator between electrodes 
W is the width of the generator 
L is the length of the generator 
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FIGURE 9. CALCULATED VELOCITY PROFILE OF TURBULENT FLOW 
IN THE MHD CHANNEL 

The velocity is the ratio of a local velocity 
based on the velocity U c at the axial center 
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FIGURE 10. CALCULATED TEMPERATURE PROFILES AT THREE DIFFERENT 
LOCATIONS FROM THE ENTRANCE OF THE MHD CHANNEL 

The constant wall temperature of 1000 k was used in 
the calculation 
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FICURE 11. CALCULATED ELECTRICAL CONDUCTIVITY PROFILES AT FOUR 
DIFFERENT LOCATIONS FROM THE ENTRANCE OF THE MHD 
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TABLE I. PLASMA MHD (CONCLUDED) 
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TABLE II. LM MHD GENERATOR 
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TABLE II. LM MHD GENERATOR (CONCLUDED) 
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TABLE III. ALKALI METAL PARAMETERS 



RESONANCE RADIA- 

IONIZATION 

PERCENT 

LENGTH FOR 90% 

GAS 

TION WAVELENGTHS 

POTENTIAL 

IONIZATION 

ABSORPTION 


(A) 

(eV) 

T=2500 K, P=760 torr 

T=2500 K, P=760 torr 

POTASSIUM 

7665 

7669 

4.45 

3.3 

52 cm 

CESIUM 

8521 

8944 

3.87 

12. 

50 cm 

SODIUM 

5890 

5896 

5.12 

0.7 

50 cm 


33 







TABLE IV. LASER-PLASMA INTERACTION 
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PARAMETERS NECESSARY FOR CALCULATING THE MHD GENERATOR 
PERFORMANCE 


NO. 

PARAMETERS 

SYMBOL 

CODE 

REMARKS 

1 

TYPE OF GAS 



NECESSARY FOR 

2 

TYPE OF SEED MATERIAL 



CALCULATING THE 

3 

COMPOSITION RATIO, SEED/GAS 



STATE VARIABLES 

A 

AMOUNT OF SEED 



AND THE ELECTRICAL 

5 

GAS DENSITY 

N 
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CONDUCTIVITY 
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SEED GAS DENSITY 
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MOLE FRACTION OF SEED 

F 



8 

INITIAL ELECTRON DENSITY 
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INITIAL ELECTRON TEMPERATURE 
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GEOMETRY 

14 

DUCT WIDTH 

x or W 



15 

DUCT HEIGHT 

y or L 



16 

MAGNET LENGTH 

b 



17 

DUCT WALL TEMPERATURE 

T 

w 


NECESSARY FOR 

18 

MAGNETIC FIELD INTENSITY 

B 


CALCULATING THE 

19 

FLOW RATE 

m 


MHD OUTPUT POWER 

20 

GENERATOR COEFFICIENT 

K 



21 

LOAD CURRENT 

I 



22 

INTERACTION COEFFICENT 
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TABLE VI. LASER-PLASMA INTERACTION PARAMETERS USED 
FOR LASER-DRIVEN MHD 


LASER 

10. 6 v C0 2 PULSED 

WORKING FLUID 

AIR - 1 atm 

BREAKDOWN 


THRESHOLD 
LASER POWER 

1.3 x 10 W/cm 

LASER POWER 

100 - 150 J, 200 ns 

10 Hz (EQIUV . CW MODE 50 - 75 MW) 

AVERAGE MEDIUM 
TEMPERATURE IN 
THE GENERATOR 
VOLUME 

> 1600 K 

ELECTRON 

DENSITY 

, n 16 _ in 19 -3 
10 10 cm 

PROPAGATING 

VELOCITY 

10^ - 10^ cm/s 

PLASMA LENGTH 

15 cm 

PLASMA RADIUS 

1 cm 

SEED MATERIAL 

CESIUM 
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TABLE VII. COMPARISON OF EFFECTS OF TURBULENCE AND 
NON-EQUILIBRIUM IONIZATION 


PARAMETERS 

TURBULENCE AND 
NON-EQUILIBRIUM 
IONIZATION 
IGNORED 

TURBULENCE AND NON-EQUILIBRIUM 
IONIZATION TAKEN INTO ACCOUNT 

geometry 

rectangle cross 
section 

rectangle cross section 

, , 3 

n , 1/cm 
e 

19 

1.315x10 

19 

1.315x10 

V K 

2500 

2500 (3100 *) 

2 

j , A/ cm 

0.3360 

0.3180 

K 

0.231 

0.2843 

B 

3.17 

2.31 

p , W/cnf* 

out 

0.225364 

0.225364 

p. . W/cnf* 
r m 

**0.420961 

0.521918 

n, % 

53.5 

43.2 


* The temperature necessary to achieve the efficiency of 53.5% 

** Based on T , n and U, P. = (n icT )AU per unit volume where 
e e in e e r 

A is the cross section area 
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